Vortex-lattice pinning in two-component Bose-Einstein condensates 
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We investigate the vortex-lattice structure for single- and two-component Bose-Einstein conden- 
sates in the presence of an optical lattice, which acts as a pinning potential for the vortices. The 
problem is considered in the mean-field quantum-Hall regime, which is reached when the rotation 
frequency fl of the condensate in a radially symmetric trap approaches the (radial) trapping fre- 
quency u> and the interactions between the atoms are weak. We determine the vortex-lattice phase 
diagram as a function of optical-lattice strength and geometry. In the limit of strong pinning the 
vortices are always pinned at the maxima of the optical-lattice potential, similar to the slow-rotation 
case. At intermediate pinning strength, however, due to the competition between interactions and 
pinning energy, a structure arises for the two-component case where the vortices are pinned on lines 
of minimal potential. 

PACS numbers: 03.75.Kk, 67.40.2w, 32.80.Pj 
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I. INTRODUCTION 

A characteristic property of a superfluid is that it 
supports angular momentum through quantized vortices 
0, Q . If many vortices are present in the system, they 
arrange themselves in a hexagonal ( Abrikosov) lattice [3| . 
The latter was observed experimentally in type-II super- 
conductors [HIB] and more recently in a single-component 
Bose-Einstein condensate (BEC) @, 0- In the case of ro- 
tating two-component condensates, the ground-state vor- 
tex lattice can have a non-hexagonal structure, since the 
vortices in both components can move with respect to 
each other. It was predicted that, among other struc- 
tures, the vortex lattice in a two-component condensate 
can have an interlaced square structure m. This struc- 
ture was indeed observed experimentally [9(. 

The Abrikosov lattice structure arises in a single- 
component BEC due to vortex-vortex interactions. One 
can also apply an optical-lattice potential to the conden- 
sate consisting of a regular pattern of potential minima 
and maxima. The prediction [l(| [TT|, [HI that vortices 
are pinned at the maxima of an optical-lattice poten- 
tial for sufficiently large strength was indeed confirmed 
experimentally |13|. Other theoretical work on rotating 
Bose gases in optical lattices focuses on the system near 
the superfluid-Mott-insulator transition (13 . fl5l . Il6j and 
on the effect of incommensurability between the vortex 
lattice and the optical lattice [171 ] . 

A remaining challenge is to determine the phase di- 
agram of a rotating two-component condensate in an 
optical-lattice potential. Previous theoretical work fo- 
cused on the slow-rotation limit and neglected the effect 
of the non-zero total particle density inside the vortex 
core [TTI ]. In the interlaced vortex-structure in a two- 
component BEC, the vortex lattices do not lie on top of 
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each other but are displaced by a fixed offset. This fea- 
ture shows that the inter-component interaction inside 
the vortex cores is important. Therefore, the approach 
from Ref. [ll[ is not suitable to determine all the vortex 
structures in rotating two-component BEC's. 

Here, we consider instead the mean-field Quantum Hall 
regime, where the angular momentum of the condensate 
is so high that the wave function resides in the low- 
est Landau level (LLL), but mean- field theory remains 
valid [H, Eil • In this regime the wave function is com- 
pletely determined by the positions of the vortices and 
no further approximations are needed. We extend the 
method from Refs. H [l8[ to calculate the optical-lattice 
energy for a wave function in the LLL. The result is used 
to determine the phase diagram of a single-component 
condensate in an optical lattice of arbitrary geometry 
and a two-component condensate in a square optical lat- 
tice. For single-component condensates, we find phases in 
which the vortices are pinned on lines of maximal poten- 
tial and phases in which they are pinned at the pinning 
centers, which is consistent with previous theoretical and 
experimental results 0, @, 0, 03, lH GJ, 0- In the two- 
component case, we find the interlaced square lattice in 
the absence of pinning, and also a new phase where the 
vortices are pinned on lines of minimal potential. 

The remainder of this paper is organized as follows. 
In Sec. [Ill we evaluate the energy functional, including 
the optical lattice, of the system in the LLL regime. The 
result is used in Sec. IIIH to determine the vortex phase di- 
agrams of single- and two-component BEC's in an optical 
lattice. In Sec. |lV]we present our conclusions. 



II. VORTEX PINNING IN THE LLL: THEORY 

We consider a rotating two-component BEC in an 
optical-lattice potential. In this section, we discuss the 
condensate wave function, the single particle energies, 
and interaction energies in the LLL regime (see also 
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Refs. |8|, [ill). Next, we evaluate the contribution of the 
optical-lattice potential to the energy of the system. 



A. Energy functional 

We are mainly interested in the two-dimensional (2D) 
ordering of the vortices. Therefore, we assume that the 
condensate has a small effective size d z in the z-direction 
and consider a 2D wave function. At the mean-field 
level, a two-component BEC is described by two macro- 
scopic condensate wave functions ip\ (r) and 1^2 (r) , where 
r = (x, y). The condensate rotates with an angular veloc- 
ity £1 around the z-axis, thus it is convenient to transform 
to a frame that co-rotates with the condensate. The wave 
functions ipj describing the condensate are found by min- 
imization of the energy functional in the rotating frame 



AC = 



E 



dv^(r)(hj-QL z )i>j(r) 



Vj+Vol, (1) 



where ipj is normalized to the number of particles of 
species j, Nj. The single-particle Hamiltonians hj are 



given by hj 



-(n 2 /2Mf)V 2 + MjwV/2, with Mj the 



mass of a particle of species j, V = (d x ,d y ), r = |r|, and 
u! is the frequency of the magnetic trapping potential in 
the radial direction, which we assume to be the same for 
particles of both species. The angular momentum oper- 
ator is L z = —ihz ■ (r x V). The interaction energy V/ 
in Eq. ([T]) reads 

Vi = j-Jdv (||^(r)| 4 + ||^ 2 (r)| 4 

+ .9i 2 |Mr)| 2 |^(r)| 2 ), (2) 



where the intra- and mier-component interaction 
strengths are given by gj — Anh 2 aj /Mj and g\2 — 
2irh 2 ai2/Mi2, respectively. Here, Mi 2 is the reduced 
mass, and aj and Q42 denote the intra- and inter- 
component scattering lengths, respectively. In this work 
we consider only repulsive interaction: (71,(72,512 > 0. 
The optical-lattice potential couples to both species in 
the same way and its energy Vol in Eq. {!]) is given by 



OL 



rfrF L(r)(|Vi(r)| 2 + |^(r)| 2 ), (3) 



where the optical-lattice potential is given by 

Vol (r) = V [ cos (k x • r) + cos (k 2 ■ r)] , (4) 

with the kj's denoting the reciprocal optical-lattice vec- 
tors, and Vo the optical-lattice strength. The relation 
between the real-space optical-lattice vectors and the 
ki's is ki = (27r/AoL)b2 X z and k 2 = -(27r/AoL)bi x z, 
where Aol = |bi X b2 1 is the area of the unit cell of the 
optical lattice. We restrict our analysis to optical lat- 
tices which have a rhombus-shaped unit cell. The angle 
between the lattice vectors is denoted by (j>. 



B. Lowest Landau level 

It is shown in Ref. [HI that in the absence of an op- 
tical lattice (Vo = 0) and in the presence of weak in- 
teraction (V/ small) the system enters the LLL regime 
when Q | u>. This result can be easily extended to the 
case where Vol is small. More specifically, the criterium 
is that interactions and the optical lattice do not cause 
transitions to higher Landau level states, i.e., it must 
hold that gjUj, (?i2v /n i n 2 ^ hco and Vo *C Hu>, where hu 
is the Landau level gap and rij is the density of species 
j particles. In this regime, the macroscopic wave func- 
tion ipj is completely determined by the position of the 
vortices and is given by 



(5) 



where lj = y/h/MjU) is the magnetic length, z = x + iy 
is the complex position coordinate, Aj is a normaliza- 
tion constant, and are the complex positions of the 
vortices in the condensate of species j. 

We assume that the vortex positions in Eq. J5| form an 
infinite regular 2D lattice, and we project the functionals 
Vj and Vol onto the space of wave functions with this 
property. The vortex lattice in condensate 1 is spanned 
by lattice vectors Cj, which are parameterized as 



.4 



ci 



VL 



p sml 



(1,0) and C2 = \/pAvi, sin #(cot 9, 1), 



where p = |c 2 |/|ci| is the ratio between the lengths, the 
angle between the Cj's is 9, and Ci lies along the x-axis. 
Since, for a given lattice, the lattice vectors can always 
be chosen such that p > 1 and 7r/3 < 6 < tt/2, we restrict 
p and 9 to these ranges. The area of the vortex-lattice 
unit cell is Ayh = |ci XC2I. The reciprocal lattice vectors 
are Ki = (27t/Avl)c2 X z and K 2 = — (27t/Avl)ci x z, 
where z is the unit vector in the z-direction. Note that 
within our approach it is not necessary that ^4ol = ^vl- 
The vortex positions in the condensate 1 are given by 
5 = {m,c, + vo}, where the m, are integers, repeated 
indices i £ {1,2} are summed over, and Vo is the offset of 
the vortex lattice in condensate 1 from the origin. Then, 
it holds for the particle density that 



|V>i(r)| 2 = /(r-v ) e 



-r /<y l 



1 1 

Wlth 72=£2- 
°1 l l 



A 



VL 



where <7 X is the effective condensate size in the radial 
direction and / is a structure function that is zero at the 
positions {micA and has the lattice periodicity: f(r + 
m iCi )=f(r) dll|. 

In this paper, we want to study the situation in which 
the vortex lattices in both components are commensu- 
rate, i.e., in which they have the same geometry and unit 
cell area. It is therefore natural to assume that atoms 
from species 1 and 2 are similar and to restrict our anal- 
ysis to the case where Ni ~ N%, Mi ~ M2, and gi ~ g 2 . 
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In the remainder, we drop the subscripts for these quan- 
tities. Thus, we assume that the vortex lattices in both 
components are the same, up to a constant offset ro: If 
the set of vortex positions in component 1 is S, the set 
of vortex positions in component 2 is 5 + ro. Then, both 
components are described by the same structure function 
/, have equal effective radial size a = <ti, and equal vor- 
tex lattice unit cell area Avl- Note that in contrast to 
the non-rotating case, the system will not phase-separate 
when gi2 > JQ\9%- The reason for this is that the system 
is already effectively phase-separated (albeit incomplete), 
when it has an interlaced vortex-lattice structure Q . The 
particle density in component 2 is given by 



|^ 2 (r)| 2 = /(r - v„ - r„)e- 



2 /a 2 



When is close to u, the density spreads out in the radial 
direction, so that a 2 /Avl ^> 1, and it can be shown Q 
that the energy functional reads 

^ 2Nh(u - tt)a 2 gN 2 , 
IC = y — — J —+ J - {I + £12/12) + Vol, (6) 



P 



2Tra 2 d? 



where 512 = 912/9- The first term in Eq. ([6]) is the 
contribution from the single particle Hamiltonians hj. 
The quantities I and -Z12 describe the intra- and inter- 
component interaction energy, respectively, and are given 
by 

/ = ^|/k| 2 and J 12 =^|/ K | 2 cos(K-r ), (7) 



K 



K 



where the /k are the Fourier coefficients of the function 
/(r) and are given by 

for a reciprocal lattice vector K = miKi + J712K2, with 
mi and 777,2 integers. In the next section we calculate 
Vol- 



C. Optical-lattice energy 



Using the Fourier expansion of the function /(r), we 
normalize the particle densities |-0j(r)| 2 to N: 

\ M r)\ 2 ^JLy2 fKe ^ vo)e -rV^ and 



\ M r)\ 2 = JL^/ Ke «-(r-vo-r„ ) e -, 2 /. 2 ; 
™ K 

where we defined / K = /k/Ek- /k'^*' ""^' 7 ^' 2 ' 4 
and / K = /K/EK'/K'e" ,;K '- (ro+Vo) e-' j2K ' 2 / 4 . We use 
that when — > w, it holds that er 2 /AvL ^ 1 and since 
K 2 ~ 1/Avl if K 7^ 0, only those terms in the summa- 
tions in the denominators of /k and /k for which K = 



survive. Hence, /k = /k — /k- By substituting the 
expressions for | , 0i(r)| 2 and |-02(r)| 2 into Eq. we find 



Vol = 



NVosr* k ■ v , 
2 ^ 



1 



-iK-r, 



where 



dv e iK r e- r 



^-(K+k,) 2 ^ 



7* 2 



cos (kj • r) 



We again use that in the fast-rotating limit (K + kj)<7 2 ^> 
1 or (K — kj)cr 2 ^> I unless K = — kj or K = kj, respec- 
tively. Thus, the Gaussian terms of Gj,K are very small 
unless their argument is zero. Therefore, it is reason- 
able to approximate GjK by the sum of two Kronecker 
delta's: Gj,K = <^K,kj +<5-K,k j - Using that Gjjz and /k 
are even in K, we find that 

Vol = NV /k^K.k, {cos (kj ■ v ) 

+ cos [kj • (v + r )]}. (9) 

From Eq. ([9]) it follows that the system can only gain 
pinning energy if there are reciprocal vortex-lattice vec- 
tors K equal to ki and/or k2. In that case the vortices 
are pinned on equally spaced lines, as we show below. 
Assume for concreteness that ki = K = miKi + 7712X2, 
with mi and 7712 integers. Using the definitions of the 
reciprocal optical-lattice and vortex-lattice vectors given 
above, this is equivalent to 777iCi—m2C2 = (AvL/^OL)b2- 
Since the Cj's are linearly independent vectors, there is a 
pair of real numbers (7*1,^) such that riCi + 7'2C2 = bi. 
By taking outer products between the left- and right- 
hand sides of the last two equalities we obtain | mi 7*2 + 
m27*i| = 1. By inverting the expressions for the bj's, 
we reach the conclusion that the coefficient of bi in the 
expansion of both Cj's is an integer. Thus, the vortices 
are pinned on the collection of lines {nhi + r\>2 + Vrj, n G 
Z and r € R}. Analogously, we can show that if there is a 
reciprocal vortex lattice vector K equal to k2 the vortices 
are pinned on the lines {rbi+nb2+vo, 77 € 1 and r € R}. 
Since vortices are density minima, one intuitively expects 
that the lines are always lines of maximal potential and 
that Vo = 0. Nonetheless, it turns out that in the two- 
component case there are regions in the phase diagram 
that have non-zero Vo- 

In the case that there are reciprocal vortex-lattice vec- 
tors K Q and Kj,, equal to ki and k 2 , respectively, the 
vortex positions lie at the intersection of the two collec- 
tions of lines we mentioned above, i.e., on the positions 
{mibi + m 2 b 2 + v , m, € Z}. 



III. VORTEX PINNING IN THE LLL: PHASE 
DIAGRAMS 

In this section we determine phase diagrams of single- 
and two-component BEC's using the results for the en- 
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FIG. 1: Phase diagram for a single component BEC in an 
arbitrary optical-lattice potential. On the horizontal axis, 
we show the angle 4> of the optical lattice, from <j> = n/3 
(hexagonal lattice) to (f> = n/2 (square lattice). On the ver- 
tical axis, we show the optical-lattice strength Vo in units of 
fj,o = gN/ira 2 dz. The unpinned, line-pinned, and fully- pinned 
phases are indicated by UP, LP, and FP, respectively. 



ergy functional in Eqs. (J5][5]) and the optical-lattice en- 
ergy Eq. ijSJ). We consider first the case of a single- 
component condensate in an optical lattice with arbitrary 
unit-cell angle and then the case of a two-component con- 
densate in a square optical lattice. For simplicity, we re- 
strict our analysis to vortex lattices with one vortex per 
optical-lattice unit-cell, i.e., to the case Avl = ^.ol = A. 
Then, the first term in Eq. ^ is a constant, which we 
drop. 



A. Single-component lattices 

For the single-component case, the energy JC S of a con- 
densate of N particles in an optical lattice is found from 
Eq. ([6]) by setting c/12 and i"o to zero and dividing all 
terms by a factor 2 



K K = 



9N 2 
4ira 2 d 



■I + NV faMMj cos ( k j ' v o) 



In Fig. [T] we display the phase diagram, which is ob- 
tained by numerical minimization of K s as a function of 
the vortex lattice parameters for given optical lattice an- 
gle 4> and strength Vo. The latter is represented in units 
of /io = gN/ira 2 d z , which is equal to the chemical poten- 
tial fj,, up to a numerical factor of order 1. The unpinned, 
line-pinned, and fully-pinned phases are indicated by UP, 
LP, and FP, respectively. We describe these phases be- 
low. 

In the unpinned phase, there is no K equal to ki or 
k2 and the last term of K. s is zero. The vortex lattice 
ignores the optical lattice and the relative orientation of 
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FIG. 2: The line-pinned and fully-pinned vortex lattices for 
the case of an arbitrary optical-lattice angle <f> and a square 
optical lattice, respectively. The dashed lines are lines of max- 
imal potential, the vortices are represented by circles. In (a) 
the pinning centers (not shown) are located on the dashed 
lines and in (b) they are represented by the crosses. 



the two lattices is not correlated. The incommensura- 
bility between the vortex- and optical-lattice causes the 
optical-lattice potential energy to average to zero. The 
optimal unpinned vortex lattice has p = 1 and 9 = tt/3. 
This corresponds to the Abrikosov structure, and the re- 
sult is consistent with experimental and theoretical re- 
sults obtained outside the LLL regime [f| 0] ■ 

The line-pinned phase for an arbitrary optical-lattice 
angle <f> is shown in Fig. [DJa), where the x-axis is in the 
direction of the dashed lines. In this phase the vortices 
(circles) are pinned on the (dashed) lines of maximal po- 
tential, specified by {rbi + ni>2, n € Z and rel}. The 
pinning centers are not shown, but they are located on 
the dashed lines. It holds that ci = t»i, Vo • y = 0, and 
the lattice vector C2 is such that the line-pinned vortex- 
lattice resembles the Abrikosov lattice structure the clos- 
est: C2 = bi/2 + (b2 • y)y. The energy does not depend 
on the x-coordinate of Vo: if we shift the vortex lattice in 
the horizontal direction its energy does not change. For 
the optical-lattice angle (f> = 7r/3, the line-pinned vortex 
lattice has precisely the Abrikosov lattice structure. 

The fully-pinned phase is shown in Fig. [2jb) for the 
square optical lattice. Here, all vortices are pinned on the 
maxima of the optical lattice, which are represented by 
crosses: Ci = bi, C2 = b2, and vo = 0. The vortex lattice 
exactly matches the geometry of the optical lattice. 



If we start in Fig. [T] in the unpinned phase in the 
case of a square optical lattice (4> — n/2) and increase 
the strength Voj there is a transition to the line-pinned 
phase at Vq = 0.007/io- When we increase Vo even fur- 
ther, we find a second transition from the line-pinned to 
the fully-pinned phase at Vq — 0.018/io- The transitions 
are caused by the fact that the system can gain pinning 
energy by transforming to a state in which the vortices 
are pinned on the lines or points of maximal potential. 
For optical-lattice angles <f> < 7r/2, one observes that a 
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smaller optical-lattice strength Vb is sufficient to trig- 
ger the transitions from the unpinned to the line-pinned 
and from the line-pinned to the fully-pinned phase. At 
optical-lattice angle <f> = 7r/3, it takes an infinitesimal Vo 
to achieve orientation locking. These observations can be 
explained by the fact that for lower optical-lattice angles 
<f> the fully-pinned and line-pinned structures resemble 
more closely the hexagonal Abrikosov lattice structure, 
which is the optimal structure in the absence of pinning. 

The orientation locking of the vortex lattice in an 
hexagonal optical lattice (0 = 7r/3) was observed ex- 
perimentally [I3I ]. Contrary to our results, the au- 
thors find that a non-zero minimum pinning strength 
is needed for orientation locking, which they suggest is 
due to long equilibration times of the system. Such non- 
equilibrium effects can not be investigated with the equi- 
librium theory presented here. The fully-pinned phase 
for a square optical lattice ((f) = it/ 2) was also found 
earlier [H], LL2|, [l3| • The structural lattice transitions in 
Fig. Q] typically occur for Vo w O.Ol^o ~ 0.01/U, which 
is in rough agreement with previous theoretical predic- 
tions [U [HJ and experimental observations 13]. Note 
that our line-pinned phase is similar to the pinned phase 
found in Ref. [10| for a one- dimensional optical lattice, 
although the geometry of the optical lattice in our sys- 
tem is two-dimensional. Furthermore, we note that the 
half-pinned phase mentioned in Refs [ToL fill ] is equivalent 
to the line-pinned phase for <p = 7r/2. 



B. Two-component lattices 

The energy functional K, of a rotating two-component 
condensate in an optical lattice is given by Eqs. (6-9). In 
Ref. Q, the system was analyzed in the absence of pin- 
ning, where Vo = 0. The authors find that the structure 
of the ground-state vortex-lattice depends on the value 
of <7i 2 . We show two possible configurations for Vb = 
in Fig. [3J where the empty (filled) circles correspond to 
vortices in component 1 (2). The interlaced square struc- 
ture, mentioned in the introduction, is the ground state 
when 0.37 < 512 < 0.93 (see Fig.GH»). For g 12 > 0.93, 
the vortex lattice deforms into a rectangular structure, as 
shown in Fig. [3^b) . The ratio between the lengths of the 
two vortex lattice vectors p continuously increases with 
increasing gi2- For further details we refer the reader to 
Ref. 1. 



The vortex phase diagram of a two-component conden- 
sate in the presence of a square optical lattice is shown in 
Fig. [U On the horizontal axis we show the scaled inter- 
component interaction strength 1712, and on the vertical 
axis the optical-lattice strength V"o in units of /^o- First, 
we notice that the left below corner of the phase dia- 
gram contains unpinned phases which are described in 
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FIG. 3: The interlaced square and rectangular vortex lattices 
in a two- component condensate in the absence of pinning. 
The empty (filled) circles correspond to vortices in component 
1 (2). For further details we refer the reader to Ref. [8|. 
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FIG. 4: Vortex phase diagram of a two-component condensate 
in the presence of a square optical lattice. On the horizontal 
axis, we show the scaled inter-component interaction strength 
§12, and on the vertical axis, the optical-lattice strength Vo 
in units of [io- 



Ref. [8( and a line-pinned phase which is a trivial exten- 
sion to the two-component case of the single-component 
line-pinned phase shown in Fig. UK a). We will not dis- 
cuss these phases further. At larger values of gi2 and/or 
Vo novel phases arise. In the pair-pinned (PP) structure 
(see Fig. |3Ja)), the vortex lattice consists of pairs of a 
component 1 and a component 2 vortex, which form a 
square lattice. The vortices in a pair are displaced by 
ro = r(ci + C2) with respect to each other and the pairs 
are pinned at the optical-lattice maxima in such a way 
that the vortices in a pair lie at equidistant positions 
from their pinning center. In Fig. [5ja), we encircled one 
of these vortex pairs for clarity. The value of r at Vo = 
is 0.5, so that the vortex lattice has the interlaced square 
structure shown in Fig. [3](a). When we increase Vb, the 
value of r decreases continuously from 0.5 to 0: The sys- 
tem gains pinning energy if the vortices in a pair move 
closer together towards the pinning center. When r = 
the system is in the fully-pinned phase (FP), where the 
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FIG. 5: The pair-pinned and line-pinned two-component vor- 
tex lattices. The crosses indicate the location of the pinning 
centers, the empty circles the vortices of in component 1, and 
the filled circles those in component 2. The horizontal dashed 
lines are lines of maximal potential. In (b) the pinning centers 
(not shown) are located on the dashed lines. 

vortices in a pair are pinned on top of each other at the 
potential maxima, analogous to the fully-pinned single- 
component vortex lattice in Fig.^b). 

At low values of Vq and for gi2 « 1.2 and g±2 ~ 3, the 
system is in the unpinned interlaced rectangular struc- 
ture (see Fig.[3Jb)), which we denote by UP. At the value 
gi2 = 2.19 = <?*, the ratio p between lengths of the lat- 
tice vectors of the unpinned rectangular lattice is equal to 
4. Then, the rectangular vortex lattice is commensurate 
with the square optical lattice and can be line-pinned 
by an arbitrarily weak optical lattice. For values of g±2 
close to g*, the transition to the line-pinned phase oc- 
curs at small, but finite values of Vq . In terms of lattice 
vectors, this phase is specified by Ci = bi/2, C2 = 2b2, 
and ro = (ci + C2)/2. Interestingly, the lines on which 
the vortices in this phase are pinned are not the lines of 
maximal potential, as one would expect intuitively, but 
the lines of minimal potential, which is opposite to the 
single-component case in Fig. [UJa). The reason is that 
the total density is not minimal on the (horizontal) lines 
on which the vortices lie, but halfway between those lines. 
It follows from Eq. ([5]) that the healing length of vortices 
in the LLL regime is equal to the distance between the 
vortices and this causes the above-mentioned location of 
the minima in the total density. The minima in the to- 
tal density are pinned on lines of maximal potential, as 
expected. 



potential in the LLL regime. We have incorporated the 
effects due to the optical lattice and determined the phase 
diagram of a single-component condensate in an optical- 
lattice potential with arbitrary unit-cell angle and of a 
two-component condensate in a square optical-lattice po- 
tential. For the single-component case, we find, among 
others, phases that are pinned on lines of maximal poten- 
tial or at the potential maxima. In the two-component 
case, we find a phase in which pairs of vortices are pinned 
at the potential maxima, and a phase where vortices 
arc pinned on the lines of minimal potential. Note that 
Ref. [l6| also points out the possibility of pinning at the 
minima of the optical-lattice potential, albeit in a differ- 
ent regime than we consider here. 



As mentioned before, the criteria for the validity of 
the LLL assumption are that the interaction energy per 
particle gn and the optical-lattice strength Vq are much 
smaller than the LLL gap huj. Full pinning typically oc- 
curs at Vq — O.Ol/io (Vq = 0.1/io) for single-component 
(two-component) condensates, so that at these typical 
values (Vq <C /io ~ 9 n ) the LLL condition remains 
valid. This difference in optical-lattice strength comes 
about because of the difference in vortex filling factors 
considered. We don't expect fundamental obstacles for 
observing vortex pinning in two-component condensate, 
since the rotation velocities for which vortex pinning in 
a single-component BEC and the interlaced square lat- 
tice structure in a two-component BEC without pinning 
were observed are similar [ 9 |, 1 131 . Our results are con- 
sistent with theoretical [Tol and experimental [l3| 
work for single-component condensates in an optical lat- 
tice (ToL [Til Fy| and two-component condensates with- 
out pinning [9(, which were performed outside the LLL 
regime. Hence, we expect that our results are, at least 
qualitatively, valid also outside the LLL regime. 



A possibility for further research is to relax the restric- 
tion that there are one or two vortices per optical-lattice 
unit-cell and study the effect of incommensurability be- 
tween the vortex lattice and the optical lattice. In partic- 
ular, it is interesting to investigate the situation in which 
the unit cell of the optical lattice is smaller than the crit- 
ical vortex-lattice unit-cell size tt£ 2 , which would require 
an extension of the formalism presented here. 



IV. DISCUSSION AND CONCLUSIONS 

In summary, we have considered vortex lattices in 
single- and two-component BEC's in a 2D optical-lattice 
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